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ABSTRACT 

We present the results of a photometric and spectroscopic survey of 321 Lyman break galaxies 
(LBGs) at z ~ 3 to investigate systematically the relationship between Lya emission and stellar 
populations. Lya equivalent widths (WL yQ ) were calculated from rest-frame UV spectroscopy and 
optical/near- inhaied/ Spitzer photometry was used in population synthesis modeling to derive the key 
properties of age, dust extinction, star formation rate (SFR), and stellar mass. We directly compare 
the stellar populations of LBGs with and without strong Lya emission, where we designate the former 
group (Wtya > 20 A) as Lya-emitters (LAEs) and the latter group (W^, ya < 20 A) as non-LAEs. 
This controlled method of comparing objects from the same UV luminosity distribution represents an 
improvement over previous studies in which the stellar populations of LBGs and narrowband-selected 
LAEs were contrasted, where the latter were often intrinsically fainter in broadband filters by an order 
of magnitude simply due to different selection criteria. Using a variety of statistical tests, we find that 
Lya equivalent width and age, SFR, and dust extinction, respectively, are significantly correlated in 
the sense that objects with strong Lya emission also tend to be older, lower in star formation rate, and 
less dusty than objects with weak Lya emission, or the line in absorption. We accordingly conclude 
that, within the LBG sample, objects with strong Lya emission represent a later stage of galaxy 
evolution in which supernovae-induced outflows have reduced the dust covering fraction. We also 
examined the hypothesis that the attenuation of Lya photons is lower than that of the continuum, as 
proposed by some, but found no evidence to support this picture. 
Subject headings: galaxies: high-redshift — galaxies: evolution — galaxies: starburst 



1. INTRODUCTION 

An increasing number of high-redshift galaxies have 
been found in the last two decades using selection tech- 
niques reliant on either color cuts around the Lyman 
limit at 912 A in the rest frame (e.g., Steidel et al. 1996, 
1999) or strong Lya line emission (e.g., Cowie & Hu 
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1998; Rhoads et al. 2000; Gawiser et al. 2006). These 
two methods, which preferentially select Lyman break 
galaxies (LBGs) and Lya-emitters (LAEs), respectively, 
have successfully isolated galaxies at redshifts up to z = 
7 (lye et al. 2006; Bouwens et al. 2008). Extensive data 
sets of LBGs and LAEs have afforded detailed studies 
of galactic clustering (e.g., Adelberger et al. 1998; Gi- 
avalisco & Dickinson 2001), the universal star formation 
history (e.g., Madau et al. 1996; Steidel et al. 1999), and 
the galaxy luminosity function (e.g., Reddy et al. 2008; 
McLure et al. 2009). While the nature of LBGs and 
LAEs have been studied at a range of redshifts (e.g., 
Shapley et al. 2001; Gawiser et al. 2006; Verma et al. 
2007; Pentericci et al. 2007; Nilsson et al. 2009b; Ouchi 



2 



et al. 2008; Finkelstein et al. 2009), the epoch around z <~ 
3 is particularly well-suited to investigation of these ob- 
jects' detailed physical properties. At this redshift, the 
prominent HI Lya line (A rost = 1216 A), present in all 
LAE spectra and a significant fraction of LBG spectra, is 
shifted into the optical, where current imaging and spec- 
troscopic instrumentation is optimized. Consequently, 
there are large existing data sets of spectroscopically- 
confirmcd z ^ 3 LBGs (e.g., Steidel et al. 2003) and 
LAEs (e.g., Lai et al. 2008), where extensive multiwave- 
lcngth surveys often complement the former and, less 
frequently, the latter. 

The mechanism responsible for LAEs' large Lya equiv- 
alent widths is not fully understood, although several 
physical pictures have been proposed (e.g., Dayal et al. 
2009; Kobayashi et al. 2010). As Lya emission is eas- 
ily quenched by dust, one explanation for LAEs is that 
they are young, chemically pristine galaxies experiencing 
their initial bursts of star formation (e.g., Hu & McMa- 
hon 1996; Nilsson et al. 2007). Conversely, LAEs have 
also been proposed to be older, more evolved galaxies 
with interstellar media in which dust is segregated to 
lie in clumps of neutral hydrogen surrounded by a ten- 
uous, ionized dust-free medium (Neufeld 1991; Hansen 
& Oh 2006; Finkelstein et al. 2009). In this picture, 
Lya photons are resonantly scattered near the surface 
of these dusty clouds and rarely encounter dust grains. 
Continuum photons, on the other hand, readily pene- 
trate through the dusty clouds and are accordingly scat- 
tered or absorbed. This scenario preferentially attenu- 
ates continuum photons and enables resonant Lya pho- 
tons to escape relatively unimpeded, producing a larger 
Lya equivalent width than expected given the underly- 
ing stellar population. To date, the distribution of dust 
in the interstellar medium has only been investigated us- 
ing relatively small samples (e.g., Verhamme et al. 2008; 
Atek et al. 2009; Finkelstein et al. 2009). 

Given the different selection techniques used to isolate 
LBGs and LAEs, understanding the relationship between 
the stellar populations of these objects has been an im- 
portant goal of extragalactic research. Recent work by 
Gawiser et al. (2006) has suggested that LAEs are less 
massive and less dusty than LBGs, prompting these au- 
thors to propose that LAEs may represent the beginning 
of an evolutionary sequence in which galaxies increase in 
mass and dust content through successive mergers and 
star formation episodes (Gawiser et al. 2007). The high 
specific star formation rate - defined as star formation 
rate (SFR) per unit mass - of LAEs (~7x 10~ 9 yr^ 1 ; 
Lai et al. 2008) relative to LBGs (~ 3 x 10~ 9 yr" 1 ; Shap- 
ley et al. 2001) illustrates that LAEs are building up stel- 
lar mass at a rate exceeding that of continuum-selected 
galaxies at z ~ 3. This rapid growth in mass is consis- 
tent with the idea that LAEs represent the beginning of 
an evolutionary sequence of galaxy formation. However, 
results from Finkelstein et al. (2009) cast doubt on this 
simple picture of LAEs as primordial objects, given that 
these authors find a range of dust extinctions (A1200 = 
0.30-4.50) in a sample of 14 LAEs at z ~ 4.5. Nils- 
son et al. (2009b) also find that z <~ 2.25 LAEs occupy 
a wide swath of color space, additional evidence that 
not all LAEs are young, dust-free objects. Furthermore, 
the assertion that LAEs are pristine galaxies undergoing 
their first burst of star formation is called into question 



by the results of Lai et al. (2008). These authors present 
a sample of 70 z <~ 3.1 LAEs, ~ 30% of which are de- 
tected in the 3.6 /xm band of the Spitzer Infrared Array 
Camera (IRAC; Fazio et al. 2004). These IRAC-detected 
LAEs arc significantly older and more massive ((£*) <~ 1.6 
Gyr, (M) ~ 9 x 10 9 ) than the IRAC-undetected sam- 
ple ({U) ~ 200 Myr, (M) ~ 3 x 10 8 M©); Lai et al. 
(2008) suggest that the IRAC-dctcctcd LAEs may there- 
fore be a lower-mass extension of the LBG population. 
Narrowband-selected LAEs are clearly marked by hetero- 
geneity, and the relationship between these objects and 
LBGs continues to motivate new studies. 

When comparing the stellar populations of LBGs and 
LAEs, it is important to take into account the selection 
biases that result from isolating these objects with broad- 
band color cuts and line flux/equivalent width require- 
ments, respectively. By virtue of selection techniques 
that rely on broadband fluxes and colors, LBGs generally 
have brighter continua than LAEs. Spectroscopic sam- 
ples of LBGs typically have an apparent magnitude limit 
of K < 25.5 (0.4L* at z ~ 3; Steidel et al. 2003) while 
LAEs have a median apparent magnitude of R ~ 27 (Ga- 
wiser et al. 2006), where 1Z and R magnitudes are compa- 
rable. Even though the majority of LBGs studied to date 
are an order of magnitude more luminous in the contin- 
uum than typical LAEs, both populations have similar 
rest-frame UV colors (Gronwall et al. 2007). Therefore, 
LAEs fainter than 1Z = 25.5 are excluded from LBG spec- 
troscopic surveys not because of their colors, but rather 
because of their continuum faintness. Given the signif- 
icant discrepancy in absolute magnitude between LBGs 
and LAEs, understanding the relationship between these 
objects can be fraught with bias. An preferable approach 
to comparing these populations is to investigate how the 
strength of Lya emission is correlated with galaxy pa- 
rameters, for a controlled sample of objects at similar 
rcdshifts drawn from the same parent UV luminosity dis- 
tribution. 

Several authors have looked at the question of the ori- 
gin of Lya emission in UV flux-limited samples (e.g., 
Shapley et al. 2001; Erb et al. 2006a; Reddy et al. 2008; 
Pentericci et al. 2007; Verma et al. 2007). Shapley et al. 
(2001) analyzed 74 LBGs at z ~ 3 and constructed 
rest-frame UV composite spectra from two samples of 
"young" {U < 35 Myr) and "old" (i* > 1 Gyr) galax- 
ies, respectively. These authors found that younger ob- 
jects exhibited weaker Lya emission than older galax- 
ies; Shapley et al. (2001) attributed the difference in 
emission strength to younger LBGs being significantly 
dustier than their more evolved counterparts. On the 
other hand, Erb et al. (2006a) examined a sample of 87 
star-forming galaxies at z ~ 2 and found that objects 
with lower stellar mass had stronger Lya emission fea- 
tures, on average, than more massive objects. In a sam- 
ple of 139 UV-selected galaxies at z ~ 2-3, Reddy et al. 
(2008) isolated 14 objects with Lya equivalent widths > 
20 A and noted no significant difference in the stellar pop- 
ulations of strong Lya-emitters relative to the rest of the 
sample. Pentericci et al. (2007) examined 47 LBGs at z ~ 
4 and found that younger galaxies generally showed Lya 
in emission while Lya in absorption was associated with 
older galaxies (in contrast to the Shapley et al. (2003) re- 
sults). Probing even earlier epochs, Verma et al. (2007) 
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Fig. 1. — Redshift distribution of the sample, where (z) = 2.99 
± 0.19. 

examined a sample of 21 LBGs at z ~ 5 and found no 
correlation between Lya equivalent width and age, stel- 
lar mass, or SFR. These authors noted, however, that 
only 6/21 of the brightest LBGs had corresponding spec- 
troscopy from which equivalent widths were estimated. 
Therefore, the lack of a correlation between Lya equiva- 
lent width and stellar populations may, in this case, have 
been masked by a small sample that was biased towards 
the brightest objects. These aforementioned investiga- 
tions have shown that there does not yet exist a clear 
picture relating stellar populations to Lya emission. 

In this paper, we present a precise, systematic inves- 
tigation of the relationship between Lya emission and 
stellar populations using our large photometric and spec- 
troscopic data set of z <~ 3 observations. As an improve- 
ment over previous studies, we approach the data analy- 
sis from multiple aspects: we compare not only the stel- 
lar population parameters derived from population syn- 
thesis modeling, but also examine the objects' best-fit 
SEDs and photometry. Furthermore, all analysis is con- 
ducted on objects drawn from the same parent sample of 
continuum-bright (TZ < 25.5) LBGs. By controlling for 
continuum magnitude, we avoid the biases of comparing 
objects with significantly different luminosities while still 
retaining the ability to comment on the nature of strong 
Lya-emitting galaxies within the LBG sample. Our con- 
clusions are applicable to both LBGs and bright (TZ < 
25.5) narrowband-selected LAEs. While we are unable 
to make inferences about the population of faint LAEs, 
our study is complete with respect to bright LAEs given 
these objects' similar colors to LBGs in the rest-frame 
UV. 

We are motivated by the following questions: how do 
the stellar populations of Lya-emitting LBGs differ from 
those of other LBGs at z ~ 3 where the Lya emission 
line is weaker (or absent altogether)? To what degree 
are galactic parameters such as dust extinction, SFR, 
age, and stellar mass correlated with Lya line strength? 
What do the relative escape fractions of Lya and contin- 
uum photons reveal about the distributions of gas and 
dust in these objects' interstellar media? 



This paper is organized as follows. In §2, we present 
details of the observations and data reduction, including 
a description of the systematic technique used to calcu- 
late Lya equivalent widths. Stellar population modeling 
is discussed in §3. The properties of objects with and 
without strong Lya emission are presented in §4 and we 
discuss how our data can be used to address several of the 
outstanding questions pertaining to the physical nature 
of LBGs and LAEs in §5. A summary and our conclu- 
sions appear in §6. We assume a standard ACDM cos- 
mology throughout with H = 70 km s _1 Mpc^ 1 , H,m = 
0.3, and £1a = 0.7. All magnitudes are based on the AB 
system (Oke & Gunn 1983) 6 , with the exception of the 
infrared passbands which are in the Vega system. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Imaging and Spectroscopy 

The data presented here are drawn from the LBG sur- 
veys of Stcidcl and collaborators, with approximately 
half of the observations described in Steidel et al. (2003, 
2004) and half from subsequent programs by the same 
authors. These surveys employed photometric preselec- 
tion in the U n GTZ passbands in a variety of fields (Reddy 
et al. 2008) to target galaxies in the redshift interval z 
~ 2-3. Follow-up optical spectroscopy of a subset of 
these galaxies, paired with supplemental near and mid- 
infrared photometry, has yielded an extensive data set 
upon which multiple studies have been based (e.g., Shap- 
ley et al. 2003; Adelberger et al. 2005a; Shapley et al. 
2005; Erb et al. 2006b; Reddy et al. 2008). 

Here, we introduce a spectroscopic and photometric 
sample of z ~ 3 LBGs. These data were photometrically 
preselected with the following standard LBG U n GlZ flux 
and color cuts: 

"K<25.5, G-K< 1.2, U n - G > G -TZ + 1 (1) 

where the U n , G, and TZ passbands sample A rcst ~ 900, 
1200, and 1700 A at z <~ 3, respectively. Object detec- 
tion, color cuts, and photometry are discussed in Steidel 
ct al. (2003). Multi-object optical spectroscopy was ob- 
tained using the Low Resolution Imaging Spectrometer 
(LRIS; Oke et al. 1995) on the Keck 1 10m telescope. The 
majority of the data (93%) were taken with the blue arm 
of LRIS (LRIS-B; McCarthy et al. 1998; Steidel et al. 
2004) , and the remainder of the data were obtained with 
LRIS prior to its blue arm upgrade in September 2000. 
The LRIS-B data were collected using 300, 400, and 600 
line mm -1 grisms, which resulted in spectral resolutions 
of A/AA = 1000, 1200, and 2000, respectively. The 400 
(600) line mm -1 grism was used for 55% (39%) of the 
observations, and the remaining <~ 6% of the LRIS-B 
spectra were obtained with the 300 line mm -1 grism. 
LRIS-B rest-frame wavelength coverage extended from ~ 
900-1500 A and a typical integration time was 3 x 1800 s. 
The data were reduced (flat-fielded, cosmic ray rejected, 
background subtracted, extracted, wavelength and flux 
calibrated, and transformed to the vacuum wavelength 
frame) using IRAF scripts. Details of the data collection 

6 AB magnitude and /„, the flux density in units of ergs s — 1 
cm -2 Hz -1 , arc related by ttiab = — 2.5 logio j v — 48.6. Con- 
versions between AB and Vega magnitudes for the near-infrared 
passbands arc as follows: K S (AB) = K s (Vega) + 1.82; J(AB) = 
J(Vega) + 0.90. 
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TABLE 1 
Spectroscopic Survey Fields 



Field Name o a <5 b Field Size N LB G C 

(J2000.0) (J2000.0) (arcmin 2 ) 



Q0100* 


01 


03 


11 


13 


16 


18 


42.9 


22 


Q0142 


01 


45 


17 


-09 


45 


09 


40.1 


20 


Q0449 


04 


52 


14 


-16 


40 


12 


32.1 


13 


Q1009* 


10 


11 


54 


29 


41 


34 


38.3 


30 


Q1217 


12 


19 


31 


49 


40 


50 


35.3 


13 


GOODS-N d *t 


12 


36 


51 


62 


13 


14 


155.3 


54 


Q1307 


13 


07 


45 


29 


12 


51 


258.7 


8 


Q1549*t 


15 


51 


52 


19 


11 


03 


37.3 


48 


Q1623*t 


16 


25 


45 


26 


47 


23 


290.0 


24 


Q1700* t 
Q2206* 


17 


01 


01 


06 


11 


58 


235.3 


39 


22 


08 


53 


-19 


44 


10 


40.5 


23 


Q2343*t 


23 


46 


05 


12 


49 


12 


212.8 


26 


Q2346 


23 


48 


23 


00 


27 


15 


280.3 


1 


TOTAL 














1698.9 


321 



a Right ascension in hours, minutes, and seconds. 
b Dcclination in degrees, arcminutes, and arcseconds. 
c Number of spectroscopically-confirmcd LBGs with rcdshifts 
z > 2.7, excluding QSOs and AGN. 
d This field is also referred to as "HDF". 
*Denotes a field with near-infrared imaging, 
t Denotes a field with mid-infrared Spitzer IRAC imaging. 



tions in at least one IRAC passband, and 34/321 objects 
(11%) have both K s and IRAC detections. 

2.2. Galaxy Systemic Redshifts 

In order to prepare the spectra for subsequent mea- 
surement and analysis, we transformed each spectrum 
into the stellar systemic frame where the galaxy's center 
of mass was at rest. To do so, we employed the procedure 
of Adelberger et al. (2003) to infer the galaxy's systemic 
rcdshift from measurements of its redshifts of both Lya 
in emission and interstellar lines in absorption. For the 
spectra that clearly exhibited a double-peaked Lya emis- 
sion feature (12/321 objects), we adopted the convention 
of setting the Lya emission redshift equal to the average 
rcdshift of the two emission peaks. This technique of 
inferring a zero- velocity center-of-mass redshift, as op- 
posed to measuring it directly, was necessary due to the 
fact that stellar lines arising from OB stars (assumed to 
be at rest with respect to the galaxy) are too weak to 
measure in individual spectra at z ~ 3. Furthermore, a 
systemic redshift could not be measured from prominent 
LBG spectral signposts (e.g., Lya or interstellar absorp- 
tion lines) as these features trace outflowing gas which is 
offset from the galaxy's center-of-mass frame by several 
hundred km s" 1 (Shapley et al. 2003). 



and reduction of both the preselection and spectroscopic 
samples are presented in Steidel et al. (2003). 

Approximately 3% of the spectroscopically-confirmed 
z <~ 3 LBGs were classified as either active galactic nu- 
clei (AGN) or quasi-stellar objects (QSOs) on the basis 
of broad lines and high-ionization emission features, re- 
spectively (Reddy et al. 2008). These objects were ex- 
cluded from the spectroscopic sample, as were galaxies at 
rcdshifts z < 2.7. The final sample, spanning 13 photo- 
metric preselection fields totaling 1700 arcmin 2 , includes 
321 objects with an average redshift of (z) = 2.99 ± 0.19 
(Table 1, Figure 1). We note that this sample is distinct 
from previous studies of z ~ 3 LBGs (e.g., Shapley et al. 
2001, 2003) in that the majority of these objects have 
corresponding near- and mid-infrared photometry. 

Near-infrared photometry in the J (A c = 1.25 /im) and 
K s (A c = 2.15 /im) bands was obtained for a subset of 
the sample (8/13 fields) using the Wide Field Infrared 
Camera (Wilson et al. 2003) on the Palomar 5m tele- 
scope. 102/321 objects (32%) were detected in K s imag- 
ing and an additional 69 objects fell on the K s images 
and were not detected. We assigned K s upper limits cor- 
responding to 3cr image depths (K s ~ 22.2 (Vega); Erb 
et al. 2006b) to these 69 galaxies. J band photometry 
was also obtained for 57/102 objects (56%) detected in 
the K s sample. Details of the data collection and reduc- 
tion of the near-infrared sample are presented in Shapley 
et al. (2005) and Erb et al. (2006b). 

Mid- infrared imaging was obtained for 5/13 fields with 
IRAC on Spitzer. Observations at 3.6, 4.5, 5.8, and 8.0 
ptm were obtained for the GOODS-N (Dickinson et al. 
2003; Giavalisco et al. 2004; Reddy et al. 2006a), Q1700 
(Shapley et al. 2005), and Q1549, Q1623, and Q2343 
(Erb et al. in preparation) fields, where 3a IRAC detec- 
tion limits ranged from 25.1-24.8 (AB). The mid-infrared 
data were reduced according to procedures described in 
Shapley et al. (2005). 112/321 objects (35%) have detec- 



2.3. Lya Equivalent Width 

HI Lya, typically the strongest feature in LBG spectra, 
is characterized by its equivalent width, W^ya, where we 
use a negative equivalent width to correspond to the fea- 
ture in absorption. We present here a systematic method 
for estimating IF Lya, taking into account the various 
Lya spectroscopic morphologies that were observed in 
the sample. In particular, this method employs a more 
robust technique than used previously to determine the 
wavelength extent over which the Lya feature should be 
integrated to extract a line flux. 

We first binned the 321 systemic- frame spectra into 
one of four categories based on the morphology of Lya: 
"emission", "absorption", "combination", and "noise". 
The spectra in the "emission" bin were clearly dominated 
by a Lya emission feature, and a small subset of this 
sample exhibited two peaks in emission. The spectra in 
the "absorption" bin were dominated by a trough around 
Lya, typically extending for tens of angstroms bluewards 
of line center. The spectra deemed to be "combina- 
tion" contained a Lya emission feature superimposed on 
a larger Lya absorption trough and the "noise" spectra 
were generally featureless around Lya, save for a possi- 
ble absorption signature whose secure identification was 
hindered by low signal-to-noise. Four example spectra, 
characterized as falling into each of these four bins, are 
shown in Figure 2. 

Each spectrum, regardless of its category classification, 
was fit with two average continuum levels, one bluewards 
(1120-1180 A; Cbiuo) and one redwards (1225-1255 A; 
Crcd) of Lya; these wavelength ranges were chosen to 
avoid the prominent Si III and Si II absorption features 
at 1206 and 1260 A, respectively. We worked with both 
the spectra and the adopted continua in f\ units (erg 
s _1 cm~ 2 A -1 ). Below, we briefly describe the procedure 
for calculating W^ya for each of the four morphological 
classification bins. 
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Fig. 2. — The Lya feature varies widely in its morphology. Four 
spectra are plotted to show representative examples of objects clas- 
sified in the "emission," "combination," "absorption," and "noise" 
bins, respectively. In order to systematically calculate Lyo equiv- 
alent width, we adopted red- and blue-side continua (horizontal 
lines from 1120-1180 A and from 1225-1255 A, respectively) and 
inferred the extent of the Lya feature (thick line below each spec- 
trum) using the methodology described in §2.3. Note that in the 
case of the "absorption" spectrum shown here, the extent of the 
Lya feature appears to extend redwards of the adopted red-side 
continuum - this difference arises because the plotted spectrum is 
unsmoothed while a smoothed spectrum was employed to calculate 
the wavelength bounds of the Lyo feature. 

Emission: 189/321 objects (59%): The wavelength of 
the maximum flux value between 1213-1221 A was cal- 
culated, as well as the wavelengths on either side of the 
maximum where the flux level intersected c ro d and Cbiuo, 
respectively. These latter two wavelengths were adopted 
as the extremes of the emission feature. In a limited num- 
ber of cases (12/188 objects), double-peaked spectra were 
individually examined to ensure that this methodology 
counted both peaks as contained within the Lya feature. 
The IRAF routine SPL0T was next used to calculate the 
enclosed flux between the two wavelength bounds. The 
enclosed flux was then divided by the level of c re( j to yield 
a measurement of WL ya in A. The level of Cbiuo was not 
used in the calculation of Wh ya due to its substantial 
diminution by the intergalactic medium (IGM). 

Absorption: 50/321 objects (16%): The boundaries of 
the Lya absorption feature were calculated in the same 



manner as those of the "emission" spectra described 
above, with the exception that the flux value between 
1213 and 1221 A was isolated as a minimum and the 
"absorption" spectra were initially smoothed with a box- 
car function of width six pixels (<~ 2.5 A) in order to 
minimize the possibility of noise spikes affecting the de- 
rived wavelength boundaries of the Lya feature. These 
smoothed spectra were only used to define the extent of 
the Lya line; the original unsmoothed spectra were used 
for the flux integration in IRAF and the enclosed flux was 
divided by c rc( j to yield W\, ya . 

Combination: 31/321 objects (10%): Objects in the 
"combination" bin were characterized by a Lya emission 
feature superposed on a larger absorption trough. The 
boundaries of the Lya feature were computed by begin- 
ning at the base of the Lya emission peak and moving 
toward larger fluxes until the smoothed spectrum (see 
above) intersected c re d and Cbiuo, respectively (the same 
technique used for the "absorption" spectra). Flux inte- 
gration and division by c rcc j were furthermore identical 
to those objects discussed above. 

Noise: 51/321 objects (16%): For these spectra domi- 
nated by noise, we adopted set values for the endpoints 
of the Lya feature based on the average boundary val- 
ues of the absorption and combination spectra - 1199.9 
and 1228.8 A. (The boundaries of the "emission" spectra 
were not included in this calculation, as the spectral mor- 
phologies of the "emission" galaxies differed greatly from 
those of the "noise" galaxies). As above, the integrated 
flux was divided by the level of c ro d to yield W^ya- 

Rest frame equivalent widths ranged from —40 A < 
W Lya < 160 A, although one object (HDF-C41) had an 
equivalent width of ~ 740 A; we attributed this outlier 
to a continuum level in the spectrum comparable with 
zero and omitted this object from further analysis. The 
median equivalent width of the sample was ~ 4 A (Fig- 
ure 3), consistent with values reported by Shapley et al. 
(2001, 2003) for z ~ 3 LBGs. 

2.4. Composite Spectra 

A composite spectrum offers the distinct advantage of 
higher signal-to-noise over individual observations. We 
accordingly constructed several composite spectra from 
our sample, using, in each case, the same basic steps 
discussed below. 

The one-dimensional, flux-calibrated, rest-frame input 
spectra of interest were stacked (mean-combined) using 
the IRAF scombine routine. Each input spectrum was 
scaled to a common mode over the wavelength range 
1250-1380 A and a small number of positive and neg- 
ative outliers (< 10% of the data) were rejected at each 
pixel position to prevent poor sky-subtraction or cosmic 
ray residuals from affecting the composite spectrum. The 
final composite spectrum was then rebinned to a disper- 
sion of 1 A pixel -1 . A composite spectrum 7 of the entire 
sample is shown in Figure 4, where several photosphcric 

7 This composite spectrum is meant to represent only the average 
of the objects in our sample, not the average of the entire z ~ 3 
LBG population. There are a variety of observational biases that 
affect that relative proportions of Lyo-cmitters and Lya-absorbcrs 
selected: large Lya emission lines contaminating the G band result 
in redder U n - G colors, scattering objects into the color section 
window (Equation 1), while Lya absorption limits the dynamic 
range in continuum magnitude over which objects are selected. We 
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Fig. 3. — The distribution of Lya rest-frame equivalent widths, 
where objects with W^Lya > 20 A (dashed line) were classified as 
LAEs. The median equivalent width of the sample is ~ 4 A, and 
an outlier at ~ 740 A (attributed to an uncertain continuum level) 
is not shown for visual clarity. 

and low- and high-ionization interstellar features are vis- 
ible in addition to Lya. 

3. STELLAR POPULATION MODELING 

As the majority of the spectroscopic sample had exten- 
sive accompanying photometric wavelength coverage, we 
conducted population synthesis modeling to derive the 
key properties of age, extinction, star formation rate, and 
stellar mass. We required at least one photometric mea- 
surement redwards of the Balmer break (A rcs t = 3646 A) 
for robust population synthesis modeling. 248/320 ob- 
jects satisfied this criterion of photometry in at least one 
near-infrared or IRAC passband, including 69 objects 
with K s upper limits (and no IRAC data). We modeled 
galaxies using Bruzual & Chariot (2003) SEDs (assum- 
ing a Salpeter (1955) initial mass function (IMF) over 
the mass range 0.1-125 Mq) and the Calzetti (2000) ex- 
tinction law derived from local starbursts, where dust ex- 
tinction, parameterized by E(B-V), was estimated from 
the latter. While the Calzetti et al. (2000) law appears 
valid on average for z ~ 3 LBGs, we discuss in §4.2 some 
caveats associated with adopting this law; we also present 
in that same section a brief discussion of our adoption of 
the Bruzual & Chariot (2003) population synthesis mod- 
els. We furthermore assumed a constant star formation 
history and solar metallicity. Recent work has suggested 
that the Salpeter (1955) IMF has too steep a slope below 
1 Mq and consequently overpredicts the mass-to-light 
ratio and stellar mass by a factor <~ 2 (Bell et al. 2003; 
Rcnzini 2006). We accordingly converted stellar masses 
and SFRs to the Chabrier (2003) IMF by dividing the 
model output values by 1.8. The modeling procedure is 
described in detail in Shapley et al. (2001); we briefly 
present a summary below. 

Firstly, for the subset of the sample where Lya fell in 
the G bandpass (2.48 < z < 3.38; 238/248 objects), we 

refer the reader to Steidel et al. (2003) for a detailed discussion of 
these biases. 



corrected the G band photometry to account for the dis- 
crepancy between the equivalent widths of the sample 
(-40 A < W Lya < 160 A) and those of the Bruzual & 
Chariot (2003) model SEDs (W Lya 10 A). Accord- 
ing to the formula in Papovich et al. (2001), we applied 
a correction in the cases where the incremental change 
in G magnitude, AG, was larger than the uncertainty on 
the original photometric measurement. This procedure 
affected 51/238 objects (21%), for which (AG) mcdian = 
0.15 magnitudes. We did not correct for possible contam- 
ination in the K s band from nebular line emission ([OIII] 
AA5007,4959, R/3 A4861), although we tested that our 
reported correlations (§4) were still robust when all ob- 
jects at z > 2.974 (such that [OIII] A5007, the strongest 
of these nebular lines, is shifted into the K s band) were 
excluded from the analysis. 

Next, for each galaxy modeled, a grid of SEDs was at- 
tenuated by dust and shifted to match the redshift of the 
galaxy. These redshifted SEDs were further attenuated 
by IGM absorption (Madau 1995) and were multiplied 
by the G, 1Z, J, K s , and four IRAC channel filter trans- 
mission curves to extract model colors. A model U n — G 
color was not calculated, as many of the objects had only 
upper limits in U n due to the significant flux diminution 
in that passband from both galactic and intergalactic HI 
absorption. These predicted colors were then compared 
with the observed colors by means of the \ 2 statistic and 
a best-fit E(B-V) and age (i*) were extracted based on 
the lowest value of x 2 - We also applied the constraint 
that t* had to be less than the age of the Universe at the 
redshift of each object. A best-fit dust-corrected SFR 
was inferred from the normalization of a galaxy forming 
stars at a rate of 1 Mq yr _1 . Stellar mass, m s tar, was de- 
fined as the integral of the SFR and the age of the galaxy. 
We did not correct stellar masses for interstellar medium 
(ISM) recycling, whereby a mass fraction of stellar ma- 
terial is returned to the ISM via winds and supernovae 
(Cole et al. 2000). Given that we are concerned with 
the relative stellar populations of objects modeled in an 
identical fashion, the constant mass factor introduced by 
assuming ISM recycling is unimportant. 

Best-fit SFR, stellar mass, E(B-V), and age values 
were extracted for the 179 galaxies without photomet- 
ric upper limits. The additional 69 galaxies undetected 
in K s imaging were modeled by adopting 3a upper limits 
in 1Z - K s color as photometric data points. We tested 
that an upper limit in 1Z — K s color was a robust proxy 
for an upper limit in both age and stellar mass by per- 
turbing the 1Z — K s colors in several increments, both 
redwards and bluewards, and re- modeling the perturbed 
SEDs (holding all other colors constant). The expected 
trend that a redder galaxy would be best fit with an 
older age and a larger stellar mass was borne out for 
the entire sample. We additionally found that E(B-V) 
and SFR were generally insensitive to perturbations in 
1Z — K s color. Therefore, given a galaxy with an upper 
limit in 1Z - K s color, its model age and stellar mass were 
adopted as upper limits and its model E(B-V) and SFR 
were assumed to be best-fit values. 

Considering the entire sample of 248 objects, we found 
median best-fit SFR, stellar mass, E(B-V), and age val- 
ues of 37 M© yr" 1 , 7.2 x 10 9 Mq, 0.170, and 320 Myr, re- 
spectively. When objects with K s upper limits (and cor- 
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responding upper limits in stellar age and mass) were re- 
moved from the analysis, the medians were 51 M yr -1 , 
8.3 x 10 9 M , 0.180, and 320 Myr, respectively. A dust 
attenuation described by E(B-V) = 0.170 corresponds 
to Ay ~ 0.7 magnitudes, using the Calzetti et al. (2000) 
starburst attenuation law with R' v = 4.05. Histograms 
of the best-fit values are shown in Figure 5. 

For each galaxy, confidence intervals of the best- fit stel- 
lar parameters were estimated with Monte Carlo simu- 
lations. First, each object's photometry was perturbed 
by an amount drawn from a Gaussian distribution de- 
scribed by the photometric error. Then, the galaxy was 
remodeled with these new "observed" colors, and this 
process was repeated 1000 times. These simulations re- 
sulted in estimates of the distributions of best-fit stellar 
parameters allowed by the uncertainties in the photom- 
etry. In some instances, these confidence intervals were 
not centered on the best-fit stellar population parame- 
ters derived from fitting the photometry. In all cases, we 
proceeded to adopt the best-fit values as representative 
of each galaxy's properties, and we furthermore assumed 
that the error on each parameter was described by the 
standard deviation of that parameter's confidence inter- 
val. We used the quantity <r x / (x) to express the ratio 
of the standard deviation of each parameter's confidence 
interval, cr x , to the mean of that parameter's confidence 
interval, (x). On average, using 3cr rejection to suppress 
the effect of outliers, we found that the median values 
of cr x / (x) for each of the four best- fit stellar population 
parameters (E(B-V), age, SFR, and stellar mass) were 



~ 0.4, 1.0, 0.8, and 0.6, respectively. 

4. STELLAR POPULATIONS & Lya LINE STRENGTH 

With Lya equivalent widths and best-fit stellar pop- 
ulation parameters in hand, we now turn to examining 
the relationship between Lya emission and stellar popu- 
lations in LBGs. The aim of this analysis is to investi- 
gate the physical nature of z <~ 3 LBGs by studying how 
objects with and without strong Lya emission differ in 
the fundamental parameters of age, stellar mass, extinc- 
tion, and star formation rate. Our full complement of 
data are used in this analysis, including rest-frame UV 
spectroscopy, broadband photometry, and best-fit stel- 
lar population parameters and SEDs. We also make 
comparisons between our LBG data and narrowband- 
selected LAEs from Nilsson et al. (2007) and Gawiser 
et al. (2007), although we caution that the relationship 
between LBGs and narrowband-selected LAEs is a field 
of extragalactic research unto itself; we refer the reader 
to §5.1 for a summary of the salient points of this topic. 

4.1. Statistical Tests 

We employed several statistical methods central to our 
investigation, including survival analysis techniques that 
were capable of analyzing data with limits ("censored 
data"). It was important to include limits in the statis- 
tical analysis as a non-negligible fraction of the sample 
(28%) was undetected in K s imaging. We used ASURV 
("Astronomy SURVival Analysis") Rev 1.2 (Isobe & 
Feigelson 1990; Lavalley et al. 1992), which implements 
the methods presented in Feigelson & Nelson (1985) 
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and Isobe et al. (1986). The bivariate correlation tests 
Kendall r and Spearman p were utilized, where, for each 
test, the degree of correlation is expressed in two vari- 
ables: tk (rsR.) and Pk (Psr)- The first variable repre- 
sents the test statistic and the second variable represents 
the probability of a null hypothesis (i.e., the probability 
that the data are uncorrelated). We also used routines 
from Numerical Recipes (Press et al. 1992), specifically 
the Kolmogorov-Smirnov (K-S) test between two data 
sets (kstwo), the Spearman's rank correlation (spear), 
and the Kendall's r correlation (kendll). The K-S test 
produces two outputs: D and P, where the former is the 
test statistic and the latter is the probability that the two 



samples are drawn from the same underlying population, 
spear and kendll are duplicates of the Spearman p and 
Kendall r tests available in ASURV; we employed the 
ASURV routines when it was necessary to include data 
with limits. 

4.2. Equivalent Width Versus Stellar Parameters 

Firstly, we investigated the correlation between Lya 
equivalent width and extinction, age, star formation rate, 
and stellar mass, respectively, using the Spearman p and 
Kendall t bivariate correlation tests. For this analysis, 
we used the entire sample with population modeling, in- 
cluding galaxies with upper limits in the K s passband 
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TABLE 2 

Correlation Coefficients 



Parameter Parameter 


Kendall r 


Spearman p 
r S R a (P S R) h 




Entire Sample (248 objects) 


W Lya 


E(B-V) 

Age 

SFR 

Stellar Mass 


6.585 (0.0000) 
3.655 (0.0003) 
6.515 (0.0000) 
0.358 (0.7205) 


-0.413 (0.0000) 
0.238 (0.0002) 

-0.410 (0.0000) 
0.011 (0.8630) 




Omitting K s non- 


-detections (179 o 


bjects) 




E(B-V) 

Age 

SFR 

Stellar Mass 


5.163 (0.0000) 
4.122 (0.0000) 
5.299 (0.0000) 
0.917 (0.3592) 


-0.385 (0.0000) 
0.319 (0.0000) 

-0.394 (0.0000) 
0.068 (0.3641) 



a Test statistic. 

b Probability of a null hypothesis (i.e., the probability that a 
correlation is not present). 



(248 objects). We note that while the correlations shown 
in Figure 6 are not without scatter, we have employed 
a variety of quantitative statistical analyses that have 
yielded a consistent picture of how stellar populations 
are correlated with Lya emission in z <~ 3 LBGs. Below, 
we introduce these results. 

E(B-V), age, and SFR all exhibited very strong corre- 
lations with W-^ya (Psr < 0.0002 in each case) 8 , whereas 
stellar mass was uncorrelated (Psr = 0.8630). These re- 
sults did not change when the 69 galaxies with K s upper 
limits (and corresponding upper limits in mass and age) 
were excluded from the analysis (Table 2). In order to 
test the robustness of these results while taking into ac- 
count the uncertainty on the best-fit parameters, we used 
the bootstrapping method. We created 100 bootstrap 
samples by randomly extracting one line from each ob- 
ject's 1000-entry confidence interval of stellar parameters 
(generated by perturbing the initial photometry accord- 
ing to its errors and then re-modeling the galaxy; §3) and 
then re-ran the Spearman p test on each of these 100 ar- 
tificial samples. A null hypothesis was consistently ruled 
out at between the 3 and 5a levels for the SFR and E(B- 
V) bootstrap samples and at the 3a level for the age sam- 
ple. In the case of the stellar mass sample, where the null 
hypothesis was ruled out at only the la level, the data 
do not support a correlation between W^Lya and stellar 
mass. These bootstrap samples show that even when 
accounting for the uncertainty in best-fit parameters, a 
strong correlation exists between equivalent width and 
age, E(B-V), and SFR, respectively, such that objects 
with strong Lya emission tend to be older, less dusty, 
and lower in star formation rate (more quiescent) than 
objects with weak or no Lya emission. 

As an independent test of the above results, we next 
used the rest-frame UV spectra to investigate how the 
morphology of the Lya feature was, qualitatively, cor- 
related with best-fit stellar population parameters. For 
this analysis, we limited our study to objects satisfying 
two criteria: 1 ) population synthesis modeling had been 

8 We note, however, that SFR and E(B— V) are not independent; 
the probability that these two parameters are uncorrelated is Psr 
= 0.0000. 



completed and 2) K s photometry consisted of detections, 
not upper limits 9 . Together, these criteria isolated 179 
objects. For each stellar population parameter, we di- 
vided the objects into three groups ( "tertiles" ) based on 
their best-fit value. A composite rest-frame UV spectrum 
was then constructed from each fertile according to the 
methodology presented in §2.4, for a total of 12 com- 
posite spectra (3 tertiles x 4 parameters; Figure 7). A 
strong correlation of equivalent width with age, E(B-V), 
and SFR was observed, such that stronger Lya emission 
was more prevalent in older, less dusty, and more qui- 
escent galaxies. At the same time, the strength of low- 
ionization interstellar absorption lines decreased with in- 
creasing Lya emission strength (Shapley et al. 2003). We 
note that these trends were still present when the spec- 
tra were median-combined to make to the composites (as 
opposed to mean-combined). 

Next, we reversed the independent and dependent vari- 
ables from the previous investigation and, this time, di- 
vided the sample into tertiles according to Lya equivalent 
width. The mean equivalent widths of each fertile were 
-15.5, 2.8, and 33.1 A, respectively. We then calculated 
the average stellar parameters of each fertile, in addition 
to the error on the mean (Figure 8). Consistent with the 
results of the bivariate correlation tests and the com- 
posite spectra divided by stellar population parameters 
discussed above, we found that as the average equivalent 
width varied from -15.5 A to 33.1 A, (i*) increased from 
440 Myr to 850 Myr, (E(B-V)) decreased from 0.24 to 
0.14, and (SFR) decreased from 300 M Q yr" 1 to 70 M Q 
yr _1 . A weak positive correlation between best-fit stel- 
lar mass and equivalent width was observed, although, 
given the errors on the average masses, we chose to adopt 
the stance that stellar mass was generally insensitive to 
changes in equivalent width. 

From the three aforementioned investigations, a con- 
sistent picture has emerged: strong Lya emission is asso- 
ciated with galaxies that are older, less dusty, and more 
quiescent than their counterparts with weaker Lya emis- 
sion or the line in absorption. We note that the trend 
of increasing Lya strength with decreasing E(B-V) is 
a well-known result (e.g., Shapley et al. 2003; Reddy 
et al. 2006b; Pentericci et al. 2007). Deharveng et al. 
(2008), however, found no apparent correlation between 
Lya equivalent width and UV color in a sample of 96 
LAEs at z ~ 0.3, although we caution that the transfor- 
mation from UV color to E(B-V) is sensitive to the ge- 
ometry of the emitting stars and absorbing dust and may 
therefore not be constant for a heterogeneous population 
such as LAEs (e.g., Witt & Gordon 2000; Granato et al. 
2000). The relationship between age and W^Lya observed 
in our data - older objects exhibiting stronger Lya emis- 
sion than younger objects - is in agreement with the re- 
sults of Shapley et al. (2003). The inverse correlation be- 
tween dust-corrected SFR and Wh ya that we observe has 
been noted by several authors as well (e.g., Ando et al. 
2004; Reddy et al. 2006b; Tapken et al. 2007), although 
Nilsson et al. (2007) report a ~ 2a direct correlation be- 
tween SFR and W^ya in a sample of 24 LAEs at z ~ 3.15. 

9 Objects with K s upper limits were excluded due to their cor- 
responding upper limits in best-fit ages and masses; such limits led 
to uncertainty as to the age and mass bin in which these objects 
should be grouped. 
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Fig. 6. — Plots indicating the correlation between Wi^ ya and each of the four best-fit stellar population parameters. Upper limits in age 
and mass are indicated by downward-facing arrows and the vertical dashed line in each panel delineates W-L ya = 20 A, the adopted lower 
bound in Lya equivalent width for LAEs. The probability of a null hypothesis (i.e., no correlation) from the Spearman p test is shown in 
each plot. Note the significant correlation between equivalent width and E(B-V), age, and SFR, respectively, in the sense that galaxies 
with strong Lya emission also tend to be less dusty, older, and more quiescently forming stars. M^Lya and stellar mass are not significantly 
correlated, as evidenced by the high probability of the null hypothesis. 



In contrast with the results presented here, however, Erb 
ct al. (2006a) reported a correlation between stellar mass 
and Lya strength in a sample of 87 star-forming galaxies 
at z ~ 2. These authors found that composite spectra 
constructed from the 30 least massive galaxies ((M) = 
5 x 10 9 Mq) and the 28 most massive galaxies ((M) 
= 7 x 10 10 Mq) in their sample differed widely in neb- 
ular line strength, with the former sample exhibiting a 
pronounced Lya emission feature and the latter sample 
showing a smaller Lya feature superimposed on a larger 
absorption trough. Erb et al. (2006a) attributed these 
results to the lower velocity dispersion of the interstellar 
medium in the less massive galaxies, where, following the 
arguments of Mas-Hesse et al. (2003), velocity dispersion 
and Lya escape fraction are anti-correlated. We there- 
fore conclude here that the observed correlations we find 



between W^ya and E(B-V), SFR, and age are supported 
by the work of others, whereas the lack of correlation be- 
tween stellar mass and H^Lya may either be a tenable 
result (and evidence of LBG evolution during the 1.1 
Gyr intervening z = 2 and z — 3) or otherwise an in- 
dication that further studies of the relationship between 
W\, ya and stellar mass at z <~ 3 are necessary. 

The correlations discussed above were derived assum- 
ing the Calzetti et al. (2000) dust attenuation law. While 
this relation has routinely been applied to samples of 
high-redshift galaxies (e.g., Erb et al. 2006b; Gawiser 
et al. 2006; Gronwall et al. 2007), recent work by Reddy 
et al. (2006b), based on Spitzer MIPS observations, has 
suggested that dust extinction in young (i* < 100 Myr), 
UV-selected galaxies at z <~ 2 is often overestimated by 
modeling a given rest-frame UV color with the Calzetti 
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Fig. 7. — Composite spectra assembled from objects ordered by E(B-V), age, SFR, and stellar mass, respectively, where objects with K s 
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and the median value in each fertile is printed next to its composite. Notice how strong Lye? emission is associated with older, less dusty, 
and more quiescent galaxies. There does not appear to be a significant correlation between Lya emission strength and best-fit stellar mass. 

et al. (2000) attenuation law. At z <~ 3, there have 
been no corresponding statistical studies, but Siana et al. 
(2008, 2009) have presented preliminary evidence for 
a similar trend based on two strongly lensed objects. 
On the other hand, for a single unlensed LBG at z ~ 
2.8, Chapman & Casey (2009) find good agreement be- 
tween the dust extinction inferred from the Calzetti et al. 
(2000) law and that estimated from the ratio of SCUBA 
850 /xm and rest-frame UV fluxes. Based on their find- 
ings, Reddy et al. (2006b) and Siana et al. (2008) ac- 
cordingly proposed that a steeper, SMC-like extinction 
relation (e.g., Prevot et al. 1984) may be more appropri- 
ate for young objects. 



With these results in mind, we assessed the uncertain- 
ties inherent in the choice of attenuation law. Limiting 
our analysis to the 179 objects without K s upper lim- 
its, we rc-modclcd the 70 galaxies with t+ < 100 Myr 
with a SMC-like extinction law (and retained the orig- 
inal fits assuming the Calzetti et al. (2000) law for the 
remaining 109 objects). In this mixed sample with some 
objects modeled with a SMC-like law and others mod- 
eled with the Calzetti et al. (2000) relation, we found no 
significant correlations between Wh ya and stellar popu- 
lations. We also investigated the effect of re-modeling all 
179 galaxies, regardless of age, with a SMC-like attenua- 
tion law. In this case, we observed correlations between 



12 



> 
I 

PQ 




0.10 



-20 





w 



Lya 



20 

(A) 



40 



CD 
< 




-20 







20 



40 



W (A) 



>» 300 




-20 




W 



Lya 



20 

(A) 



40 



© 

o 
O 



(0 




-20 




W 



Lya 



20 

(A) 



40 



Fig. 8. — Best-fit stellar population parameters versus rest-frame Lya equivalent width. Omitting galaxies with K s upper limits, the 
remaining 179 objects with stellar population modeling were ordered by equivalent width and divided into three groups. The average stellar 
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group). With increasing Lya strength, galaxies are, on average, older, less dusty, and more quiescent. A weak trend of increasing stellar 
mass with increasing Lya strength is observed, although this is likely not a significant correlation. 



W^Lya and E(B-V) and between M^Lya and SFR. Given 
the uncertainties in adopting the SMC law for part or 
all of our sample, another approach consists of restrict- 
ing our analysis to objects for which the validity of the 
Calzetti et al. (2000) law has not been questioned. When 
we omitted objects younger than t+ = 100 Myr from our 
original analysis assuming the Calzetti et al. (2000) dust 
attenuation law, the trends between Wh ya and E(B-V), 
age, and SFR, respectively, were still present (but only 
at a 2cr level of significance). 

Given the lack of systematic studies of the dust attenu- 
ation law in z ~ 3 objects, and the uncertainty regarding 
the age at which a galaxy may transition from being fit 
with a SMC-like law to the Calzetti et al. (2000) law, 
we choose to present all analysis on results derived from 
assuming the Calzetti et al. (2000) attenuation law. We 
acknowledge that the choice of an extinction law is a sys- 
tematic uncertainty in our analysis, yet the persistence of 



trends between Lya and stellar populations in the sam- 
ple limited to objects with t* > 100 Myr supports the 
analysis with the Calzetti law applied to the full sample. 

We now turn to discussing our adoption of the Bruzual 
& Chariot (2003) population synthesis models. While 
other models more fully take into account the Thermally- 
Pulsing Asymptotic Giant Branch (TP-AGB) stellar 
phase (e.g., Maraston 2005), a study of young star- 
forming "BzK" galaxies at z <~ 2 by Daddi et al. (2007) 
found that the Maraston (2005) and Bruzual & Chariot 
(2003) models yielded consistent M/L ratios; the LBGs 
in our sample are as least as dominated by the young 
stellar component as the Daddi et al. (2007) BzK ob- 
jects. Furthermore, a significant fraction (~ 40%) of our 
sample does not have IRAC coverage, and therefore our 
SED fitting only includes photometry out to the rest- 
frame V band (observed K s band) where the effects of 
TP-AGB stars on the derived parameters are less pro- 
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Fig. 9. — Histograms of best-fit dust attenuations, ages, star formation rates, and stellar masses for the 34 LAEs (blue line) and 145 
non-LAEs (red line), non-inclusive of objects with K s upper limits. The ordinate is the fractional percentage of the sample and the median 
value of each sample is indicated. The LAE sample is characterized by lower dust attenuation levels than the non-LAE sample; all objects 
with E(B-V) > 0.40 are non-LAEs. The median age (763 Myr) of LAEs is significantly older than that of the non-LAEs (181 Myr). 
Furthermore, the LAEs lack the conspicuous overdensity of objects in the youngest (t* < 100 Myr) age bin that characterizes the non-LAE 
distribution. While the fraction of LAEs and non-LAEs with star formation rates between ~ 80-120 A/q yr -1 is roughly the same, LAEs 
show a larger fraction of objects with small SFRs and non-LAEs exhibit a tail of high star formation rate objects. LAEs are marginally 
more massive than non-LAEs, on average, although the distributions are similar. 

Finally, we note that our population synthesis results are 
supported by empirical color differences (§4.4). In this 
sense, the choice of synthesis code is rendered less im- 
portant - there are clearly intrinsic differences between 
objects. 

In the next section, we divide the LBGs into two groups 
according to Lya strength and compare the stellar pop- 
ulations of these subgroups in order to more generally 
comment on the differences between objects with and 



nounced. Additionally, and perhaps most importantly, 
we are concerned with trends in the data (rather than 
specific stellar population values) . It is these trends that 
are preserved regardless of model choice; Lai et al. (2008) 
find that even though the Maraston (2005) models yield 
lower stellar masses and younger ages than Bruzual & 
Chariot (2003) models for two samples of objects, the 
distinction between the two samples is significant irre- 
spective of which population synthesis model is adopted. 
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without strong Lya emission. 



4.3. The Distinct Properties of Strong Lya Emitters in 
the LBG Sample 

From the perspective of examining Lya emission, the 
simplest division of our data is according to Lya equiv- 
alent width. We turn here to comparing the properties 
of objects with strong (Wh ya > 20 A) Lya emission and 
with weak Lya emission, or the line in absorption (Wh ya 
< 20 A). Given our large sample size, these two groups 
each have sufficient membership to ensure robust statis- 
tics even when we restrict our analysis to only those 
objects with both stellar population modeling and no 
K s upper limits (34 and 145 objects, respectively). We 
call these two groups LAEs and non-LAEs, respectively, 
where we differentiate these objects from LAEs isolated 
with narrowband filters by explicitly referring to the lat- 
ter as narrowband-selected objects. The motivation for 
choosing 20 A as the delineating equivalent width stems 
from the adoption of this value as a typical limit in nar- 
rowband LAE studies (e.g., Gawiser et al. 2007; Nilsson 
et al. 2009b; Pentericci et al. 2009). In order to inves- 
tigate how stellar populations vary between LAEs and 
non-LAEs, we examined the age, E(B-V), stellar mass, 
and SFR distributions of both samples (Figure 9). We 
discuss below the striking differences between the prop- 
erties of UV continuum-bright LAEs and non-LAEs. 

There is a remarkable dissimilarity between the age 
histograms of LAEs and non-LAEs: LAEs are generally 
fit with older best-fit ages than non-LAEs: LAEs have a 
median age more than four times that of the non-LAEs 
(763 Myr versus 181 Myr). LAEs furthermore lack the 
conspicuous overdensity of young objects (i* < 100 Myr) 
that characterizes the non-LAE sample. The disparity in 
ages between LAEs and non-LAEs is reflected in the K- 
S test probability of ~ 2% that the two age samples are 
drawn from the same underlying population. 

The distributions of dust attenuations for the LAE 
and non-LAE samples differ significantly, with a median 
E(B-V) of 0.13 for the LAEs and 0.20 for the non-LAEs. 
While both populations have roughly the same distribu- 
tion of E(B-V) values below 0.30, the non-LAE distri- 
bution is characterized by a tail of high extinction val- 
ues. No LAEs have E(B-V) > 0.40, whereas 12/145 
(8%) of non-LAEs do. The K-S test probability that the 
LAE and non-LAE extinction distributions derive from 
the same parent population is ~ 2%. 

The star formation rates of LAEs and non-LAEs are 
markedly different, as the LAEs have a median rate of 
25 M yr _1 and the non-LAEs are characterized by a 
median rate of 59 Mq yr _1 . The relative quiescence of 
strong Lya-emitters is in contrast to the high SFR tail 
observed in non-LAEs, where star formation rates in ex- 
cess of 500 Mq yr _1 are recorded for 22 objects 10 . There 
is a ~ 1% probability of these LAE and non-LAE star 
formation rate distributions are drawn from the same 
underlying population. 

10 15/22 of these objects, however, are younger than 10 Myr. 
While these objects have unphysically young ages (given that the 
dynamical timescale of LBGs is > 10 Myr; Giavalisco et al. 1996; 
Pettini et al. 1998, 2001), re-modeling these galaxies with the con- 
straint that t* > 10 Myr still results in best-fit star formation rates 
of several hundred Mq yr -1 . 
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Fig. 10. — Top panel: Composite rest-frame UV to near-infrared 
SEDs of the LAE and non-LAE samples, assembled from best-fit 
SEDs from population synthesis modeling (omitting galaxies with 
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and higher dust attenuations derived for LBGs in the literature. 
Bottom panel: Un-normalized spectra; note the significant flux 
discrepancy (~ 2.3 magnitudes) between the color-selected LBGs 
and the narrowband- selected LAEs. 



Unlike age, E(B-V), and SFR, the stellar mass distri- 
butions of LAEs and non-LAEs are not strongly dissim- 
ilar. The median value of the LAEs (1.1 x 10 10 Mq) is 
only 40% larger than that of the non-LAEs (7.9 x 10 9 
Mq) and a K-S test predicts that the two distributions 
have a ~ 27% chance of being drawn from the same par- 
ent distribution. 

We have shown that three of the best-fit stellar pop- 
ulation parameters - age, E(B-V), and SFR - have 
markedly different distributions in the LAE and non- 
LAE samples. LAEs are typically older, less dusty, and 
less vigorously forming stars than non-LAEs. We tested 
the robustness of these results by constructing bootstrap 
samples from the confidence intervals of best-fit stellar 
parameters (§4.2). The mean and the error on the mean 
(a / y/n, where n is the sample size) were calculated for 
each best-fit stellar population parameter in the LAE 
and non-LAE samples. This process was repeated for 
100 bootstrap samples. For E(B-V), age, and SFR, we 
found that the means of the LAE and non-LAE sam- 
ples were consistently offset by at least their errors. In 
terms of stellar mass, the LAEs samples had consistently 
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TABLE 3 
Average Photometry 



Sample ((z)) 


(K) (a) 


(G> a (a) 


(GlGM corr.) (f) 


(Ks) h 


LAEs (2.99) 


24.31 (0.02) 


24.85 (0.03) 


24.67 (0.03) 


21.24 


non-LAEs (2.97) 


24.34 (0.01) 


24.99 (0.01) 


24.81 (0.01) 


21.36 



Galaxies with K s upper limits were omitted from this analysis, as were 
objects lacking stellar population modeling. 
a G magnitudes have been corrected to account for contributions from Lya. 
b Photometry calculated from best-fit SED, for objects lacking K s imaging. 



larger values; we note, however, that the masses of the 
LAE and non-LAE samples were comparable within their 
errors. These results, in agreement with the histograms 
discussed above, support the observed correlation that 
galaxies with strong Lya emission tend to be older, less 
dusty, and more quiescent than objects with weak Lya 
emission, or the line in absorption. 

4.4. Spectral Energy Distributions 

In order to place the results of the previous section on a 
more empirical footing, we transition here from focusing 
on stellar parameters and instead examine more broadly 
the relative colors of objects with and without strong 
Lya emission. We also contrast our data with SEDs of 
narrowband-selected LAEs from Gawiser et al. (2007) 
and Nilsson et al. (2007) with the aim of investigating 
the color differences, if any, between objects isolated with 
color cuts and narrowband filters. Our analysis is focused 
on the 179 objects in the LBG sample satisfying the cri- 
teria of population synthesis modeling and no K s upper 
limits. This empirical approach is complementary to the 
stellar population modeling discussed previously, and is 
furthermore a non-parametric probe of the differences 
between objects with and without strong Lya emission. 
Such a study is necessary given the well-known degenera- 
cies among physical properties (e.g., dust extinction and 
age) derived from stellar population modeling (Shapley 
et al. 2003). 

We find stark differences between the colors of LAEs 
and non-LAEs (Table 3), where all G magnitudes have 
been corrected for the contribution of Lya (§3). Strong 
Lya-emitters have bluer G - 1Z colors than non-LAEs 
({G - n) LA E = 0.54; (G - 7?.) non — lae = 0.65) and this 
trend is consistent when a correction for IGM absorption 
(Madau 1995) is applied to the objects' G magnitudes 

((GlGM corr. — "^)lAE = 0.36; (GlGM corr. — ^)non-LAE = 

0.47). Given that non-LAEs have younger best-fit stellar 
ages than LAEs (§4.2), we expect that non-LAEs should 
have intrinsically bluer stellar continua than objects with 
strong Lya emission. The fact that we instead observe 
redder G-1Z colors in these objects is evidence that non- 
LAEs are more strongly attenuated by dust than LAEs. 
This result is consistent with the larger E(B-V) values 
derived for the non-LAEs than the LAEs (§4.2). 

LAEs and non-LAEs also differ in their 1Z - K s col- 
ors, where K s magnitudes were inferred from best-fit 
SEDs for the 78 objects lacking near-infrared photom- 
etry. We find that LAEs are characterized by (1Z — K s ) 
= 3.07, whereas non-LAEs are bluer by 0.1 magnitudes 
with (1Z — K s ) =2.98. Given that the 1Z and K 8 pass- 
bands bracket the Balmer break for objects at z <~ 3, 



a redder 1Z - K s color is indicative of a larger Balmer 
break. A more pronounced Balmer break, in turn, is 
correlated with an older stellar population (e.g., Bruzual 
& Chariot 2003). Therefore, the result that LAEs have 
redder 1Z - K s colors than non-LAEs is empirical, non- 
parametric evidence that strong Lya-emitters are older 
on average. We have accordingly shown that the differ- 
ent colors of LAEs and non-LAEs are reflective of the 
trends discerned using stellar population modeling. 

Composite best-fit SEDs of the LAE and non-LAE 
samples were assembled by directly summing the indi- 
vidual best-fit SEDs over the rest frame range 912 A-1.0 
/*m and normalizing by the number of constituent spec- 
tra (34 and 145, respectively). In Figure 10, we compare 
the composite LAE and non-LAE spectra with best-fit 
average SEDs from two studies of narrowband-selected 
LAEs. Data from Nilsson et al. (2007) and Gawiser et al. 
(2007), where the former include LAEs at z <~ 3.1 in 
the GOODS South Field and the latter encompass LAEs 
at z ~ 3.1 in the Extended Chandra Deep Field South, 
are shown both normalized and un-normalized with our 
LAE and non-LAE composite spectra. It is immediately 
apparent that the LBG-selected LAEs and non-LAEs 
are significantly brighter than the narrowband-selected 
LAEs (~ x2.3 magnitudes in the 1Z band). The LBG- 
selected objects also have redder G - 1Z and redder 7Z - 
K s colors than both the Nilsson et al. (2007) and Gawiser 
et al. (2007) data. These differences in color are consis- 
tent with LBG-selected objects being both dustier and 
older than objects isolated with narrowband filters, as 
has been suggested by some authors (e.g., Gawiser et al. 
2007; Nilsson et al. 2009b). The relationship between 
LBGs and narrowband-selected LAEs is far from clear, 
though; some authors hypothesize a continuity between 
at least portions of the two populations based on similari- 
ties in stellar mass, color, and clustering (e.g., Adelberger 
et al. 2005b; Gronwall et al. 2007; Gawiser et al. 2007; Lai 
et al. 2008; Verhamme et al. 2008) whereas others main- 
tain that the large luminosity and W^ya discrepancies 
between LBGs and objects isolated with narrowband fil- 
ters dictate separate evolutionary paths (e.g., Malhotra 
& Rhoads 2002). 

5. DISCUSSION 

Our large data set of both spectroscopically- 
determined Lya equivalent widths and stellar popula- 
tion parameters constitutes a unique sample with which 
to investigate some of the current trends reported for 
continuum- and narrowband-selected objects. We turn 
here to examining several of the outstanding questions 
pertaining to LBGs and narrowband-isolated LAEs, with 
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Fig. 11. — Absolute UV magnitude versus rest-frame W-^ ya for 1Z < 25.5, H^Lyo, > 20 A LBGs (diamonds and triangles) and the 
Gronwall et al. (2007) sample of z = 3.1 LAEs (crosses). Note that the LBG data include objects from Shapley et al. (2003) and therefore 
span a larger range of W-L ya than that presented in this paper. A strong correlation between Mtjv an< i W^Lya is observed in the Gronwall 
et al. (2007) sample while weaker correlations are seen in the sample of z ~ 3 LBGs. We attribute this discrepancy to the larger range of 
Mtjv probed by LAEs. Object Q1009-C9 is both anomalously bright (Mtjv = -22.0) and the strongest Lya emitter in our sample (Wj^ ya 
= 158.5 A), yet this galaxy does not exhibit any spectral signatures indicative of an AGN. 



the constant aim of physically motivating our work. We 
begin below with a discussion of the necessity of com- 
paring objects drawn from the same parent luminosity 
distribution. 

5.1. A Caveat: Differing Rest-Frame Luminosities 

By virtue of broadband selection techniques, LBGs 
have characteristically brighter optical continua than 
narrowband-selected LAEs (1Z < 25.5 versus R <~ 27). 
As it is critical to disentangle how luminosity and stellar 
populations arc related before attributing galaxy differ- 
ences purely to the same mechanisms that modulate Lya 
emission, it is of interest to investigate how, if at all, the 
derived properties of narrowband-selected LAEs are a 
natural result of preferentially isolating continuum-faint, 
line-bright objects. We compare our data with the work 
of several authors to explore the dependencies of Lya 
emission, stellar populations, and spatial clustering on 
broadband luminosity and we discuss the applications 
of such investigations to both LBGs and narrowband- 
selected LAEs. 

Luminosity and Equivalent Width: Recent work has 
suggested a correlation between luminosity and Lya 
equivalent width, in the sense that few luminous LBGs 
have been observed with large equivalent widths (e.g., 
Shapley et al. 2003; Erb et al. 2006a; Ando et al. 2006; 
Ouchi et al. 2008; Vanzella et al. 2009; Pentericci et al. 
2009). Several physical pictures have been suggested to 
explain this correlation, which has been noted in sam- 
ples at redshifts z <~ 3-6. Verhammc et al. (2008) and 
Vanzella et al. (2009) hypothesized that different dust 



attenuations (from a luminosity-dependent chemical evo- 
lution history) could be responsible while Ando et al. 
(2006) suggested that an enhanced presence of HI gas 
around luminous LBGs could explain the lack of large 
Lya equivalent widths in bright galaxies (consistent with 
the prediction of more massive galaxies residing in larger 
dark matter halos that are presumably richer in HI gas) . 
Some authors, however, have conversely reported that 
the correlation between luminosity and Lya equivalent 
width is tenuous, at best (e.g., Steidel et al. 2000; Verma 
et al. 2007; Stanway et al. 2007; Nilsson et al. 2009a). 
Given the ongoing debate about such a correlation, we 
examined our extensive data set to test for a deficiency 
of luminous objects with large Lya equivalent widths. 

We restricted our analysis to the 64 objects with Wh ya 
> 20 A in order to ensure completeness as a function of 
7Z magnitude, given that spectroscopic follow-up is bi- 
ased toward objects with large, positive Lya equivalent 
widths. These objects have W-^ya < 160 A and Mtjv 
ranging from -20.0 to -22.4 (neglecting corrections for 
dust attenuation). This parameter space is comparable 
to that of Ando et al. (2006), where these authors suggest 
a threshold luminosity of M1400 = -21.5 to -21.0 above 
which objects are deficient in large W^ya values. We 
find that Mtjv and IF Lya are not correlated in our sam- 
ple, where the Spearman p test statistic, tsr, is -0.104 
and the probability of a null hypothesis, Psr, is 0.4076. 
Furthermore, we see no evidence for any such threshold 
luminosity in our data; when the sample is divided into 
bright (Mtjv < -21.0) and faint (M uv > -21.0) groups, 
the median Wh ya of each group are identical (<~ 32 A). 
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These results are contradictory to those presented by 
Shaplcy ct al. (2003) who also examined a sample of z <~ 
3 LBGs. The Shapley et al. (2003) sample, when limited 
to objects with Wi, ya > 20 A, consisted of objects with 
absolute UV luminosities between -19.5 and -22.4 (ne- 
glecting corrections for dust attenuation). In those data, 
there is a significant correlation between Muv and W^ ya 
(r SR = 0.274; P S R = 0.0002). Dividing this sample into 
bright and faint subgroups, where the boundary luminos- 
ity is Mtjv = -21.0, also yielded a pronounced difference 
in the median Wh ya values of each group: My v (bright) 
= 35 A and M uv (faint) = 44 A. Given that both our 
present sample and the Shapley et al. (2003) sample are 
composed of z ~ 3 LBGs, we hypothesize that one of 
reasons our data may not exhibit a correlation between 
Muv and W-\^ ya is that our sample includes only a lim- 
ited number of high equivalent width objects. Nine ob- 
jects have W^Lya > 70 A and just one has W^Lya > 100 
A in our data set whereas, in the Shapley et al. (2003) 
sample, 45 objects have W\, ya > 70 A and 20 have W\, ya 
> 100 A. After testing that our systematic approach to 
calculating Lya equivalent width returned values consis- 
tent with those reported by Shapley et al. (2003), we 
combined our data with the 189 objects 11 from Shapley 
et al. (2003) with Wh ya > 20 A. This combined sample 
of 243 objects, with equivalent widths measured in a uni- 
form manner, was then tested for a correlation between 
Muv and W^ ya . We found only a moderate correlation 
(r SR = 0.149; P S R = 0.0202), where we hypothesize that 
the trend between Muv and W^, ya is simply not strong 
enough to be unequivocally detected in all samples. Fur- 
thermore, the correlation between Muv and W\, ya in 
these data may be masked by the relatively small range 
of Muv probed by LBGs. Narrowband-selected LAEs 
exhibit a wider range of absolute luminosities (23.0 < 
R < 27.0) than LBGs do (23.0 < K < 25.5); any trend 
between Muv and W-^ ya will consequently be more pro- 
nounced in a sample of LAEs. In order to test this state- 
ment, we analyzed a sample of LAEs from Gronwall et al. 
(2007). We divided the data into two samples, where the 
first subsample included only objects brighter than the 
LBG spectroscopic limit and the second subsample was 
inclusive of all the data. The first subsample showed 
no correlation between Muv and W\, ya (rsR = -0.015; 
Psr — 0.9185) whereas the second subsample exhibited 
a strong correlation (rsR = 0.585; Psr = 0.0000). These 
results support our hypothesis that the trend between 
Muv and W^Lya at z ~ 3 is apparent only when consid- 
ering data encompassing a large range in absolute lumi- 
nosity. In Figure 11, we present a plot of Muv versus 
W hya for K < 25.5, W hya > 20 A LBGs and the Gron- 
wall et al. (2007) sample of narrowband-selected LAEs. 
Even if one assumes incompleteness in the Gronwall et al. 
(2007) sample at the faint end, the lack of large equiv- 
alent widths among bright objects is striking (but see 
Nilsson et al. 2009a). Ando et al. (2006) report a defi- 
ciency of bright objects with large equivalent widths in 
their sample of z = 5 - 6 objects; that fact that we ob- 

11 We limit our analysis to the 179/189 galaxies satisfying the 
criterion WL ya > 20 A, measured according to our methodology 
described in §2.3; 10 objects classified as LAEs in the Shapley 
ct al. (2003) sample do not have equivalent widths above 20 A 
when measured using our systematic technique. 



serve several such objects here while those authors do 
not can be attributed to the small sample size of Ando 
et al. (2006). 

Luminosity and Stellar Populations: Understanding 
how luminosity and stellar populations are correlated 
is critical to commenting on the nature of LBGs and 
LAEs, given that the latter are commonly an order of 
magnitude fainter in optical broadband filters. Lai ct al. 
(2008) examined a sample of 70 z ~ 3.1 LAEs, ~ 30% of 
which were detected in the Spitzer IRAC 3.6 ^m band. 
These authors found that the IRAC-detected LAEs were 
significantly older and more massive ((t*) ~ 1.6 Gyr, 
(M) - 9 x 10 9 M Q ) than the IRAC- undetected sam- 
ple ((U) ~ 200 Myr, (M) - 3 x 10 8 M Q ). In addition 
to having redder colors, the IRAC-detected sample was 
also approximately three times brighter in the R band. 
This correlation between luminosity and redness, age, 
and mass prompted Lai et al. (2008) to suggest that the 
IRAC-detected LAEs may be a low-mass extension of 
the LBG population. The relationship between lumi- 
nosity and average dust attenuation was examined by 
Reddy et al. (2008), where these authors studied a sam- 
ple of UV-selected objects at 1.9 < z < 3.4 and found no 
correlation between Muv and E(B-V) over the absolute 
magnitude range 22.0 < 1Z < 25.5. While bright objects 
do not appear to exhibit a correlation between dust at- 
tenuation and luminosity, Reddy ct al. (2008) postulate 
that objects fainter than the LBG spectroscopic limit of 
1Z = 25.5 may have lower average dust attenuations than 
brighter objects. This has been confirmed by Bouwens 
et al. (2009), who examined a sample of objects at z <~ 
2-6 over a larger dynamic range (0.1L* =3 to 2L* =3 ) than 
had been previously studied. These authors reported a 
correlation between UV luminosity and dust attenuation 
such that fainter objects are bluer. Given the correla- 
tion we find between E(B-V) and Wh ya , the evolution 
of dust attenuation with UV luminosity may have im- 
plications for the fraction of LAEs as a function of UV 
luminosity (Reddy & Steidel 2009). These results have 
highlighted that examining objects over only a limited 
range of absolute magnitudes may mask an underlying 
correlation. 

Luminosity and Clustering: Several authors have 
noted the different clustering properties of LBGs and 
LAEs, where the former exhibit a mean halo occupa- 
tion of ~ 100% (e.g., Conroy et al. 2008) in dark matter 
halos with minimum masses of <~ 10 1 hr 1 M© (e.g., 
Adelberger et al. 2005b; Conroy et al. 2008) and the 
latter cluster more weakly and appear to occupy only 
1/100 - 1/10 of similarly-clustered dark matter halos 
with lower limit masses of ~ 10 10 6 h~ x M Q (e.g., Ga- 
wiser et al. 2007). Given the direct dependence of clus- 
tering strength on UV luminosity (e.g., Ouchi et al. 2004; 
Adelberger et al. 2005b; Lee et al. 2006), the spatial dif- 
ferences between LBGs and LAEs may simply reflect the 
discrepancy in typical UV luminosity between the two 
samples. It is therefore of interest to investigate how the 
clustering strength of LAEs correlates with Lya luminos- 
ity or Lya equivalent width, given that these quantities 
govern the selection of LAEs. Discerning the relative 
clustering properties of objects with and without strong 
Lya emission within the same UV luminosity range is 
critical to understanding the correlation between Lya 
and galactic properties. 
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These results - that objects more luminous in the 
rest-frame UV may 1 ) be deficient in large Lya equiv- 
alent widths and 2) be generally redder, older, and 
more massive relative to less luminous objects - may be 
the driving factors responsible for (faint, high equiva- 
lent width) narrowband-selected LAEs being typically 
younger, bluer, and less massive than (bright, lower 
equivalent width) LBGs. 

As the LBG spectroscopic limit of 1Z < 25.5 necessar- 
ily limits our sample to a smaller absolute magnitude 
range than than probed by narrowband-selected LAEs 
(23.0 < R < 27.0), the luminosity-dependent trends re- 
ported above for Lya equivalent width, stellar popula- 
tions, and clustering (e.g., Ando et al. 2006; Lai et al. 
2008; Adelberger et al. 2005b) may only become appar- 
ent in narrowband-selected samples where the absolute 
luminosity range is significantly larger than it is in sam- 
ples of LBGs. We remind the reader that while our con- 
clusions are applicable to both LBGs and bright (1Z < 
25.5) narrowband-selected LAEs, we are unable to make 
inferences about the population of faint LAEs. We refer 
the reader to Lai et al. (2008) for a discussion of these 
objects. 

5.2. LBG and LAE Equivalent Width Distributions 

The equivalent width of the Lya feature is a common 
benchmark that can be used to make comparisons be- 
tween studies of LBGs and LAEs. Observing programs 
of flux-selected Lya-emitters typically employ a narrow- 
band (~ 50-100 A) filter centered on the redshifted Lya 
line, paired with a broadband filter used to characterize 
the local continuum. Samples are defined according to 
a minimum Lya equivalent width - typically <~ 20 A in 
the rest frame - where Wi,y a is calculated photometri- 
cally by ratioing the flux in the narrowband filter to the 
flux density in the broadband filter. LBGs, on the other 
hand, are selected according to broadband flux and color 
cuts, with no implicit requirement on Why a- Equiva- 
lent widths of LBGs are most commonly calculated spec- 
troscopically, as these objects are characterized by rela- 
tively bright continua (7Z < 25.5). Given these different 
selection techniques and methods of calculating W\, ya , 
it is of interest to compare the equivalent width distri- 
butions of LBG-selected LAEs and narrowband-selected 
LAEs. A difference in the W^Lya distributions of these 
two populations would be indicative of an underlying dis- 
similarity between LBG-selected LAEs and narrowband- 
selected LAEs, even though both classes of objects can 
be broadly classified as strong Lya-emitters 12 . 

For this analysis, we combined our present sample of 
64 LBG-selected LAEs with the 179 objects from Shap- 
ley et al. (2003) satisfying the criterion of W\, ya > 20 
(measured using our methodology described in §2.3). In- 
tegrating these two data sets yields a large sample of 
2^3 objects for which equivalent widths were mea- 
sured in an identical manner. We find that these data 
are characterized by a median equivalent width of 42 A 

12 Reddy et al. (2008) found that the standard LBG color selec- 
tion criteria (Equation 1) do little to bias the intrinsic distribution 
of W^yc, for z ~ 3 LBGs, given the low probability of an object 
having a strong enough Lya feature to perturb its colors out of the 
LBG selection window. We therefore attribute any differences in 
the H^Lya distributions of LBGs and narrowband-selected LAEs 
to intrinsic dissimilarities between the two populations. 



and mean of 57 A. In contrast, a sample of 160 LAEs 
at z = 3.1 in the Extended Chandra Deep Field South 
observed by Gronwall et al. (2007) had a median of <~ 60 
A and a mean of ~ 80 A. We used the K-S test to quan- 
tify the likelihood that the LBG-selected LAEs and the 
narrowband-selected LAEs derived from the same parent 
equivalent width distribution. When the entire sample 
of objects from Gronwall et al. (2007) was considered, 
we found a probability of 3.5 x 10~ 4 that the objects 
shared a common equivalent width distribution. Alter- 
natively, when we limited the Gronwall et al. (2007) data 
to include only objects satisfying the LBG spectroscopic 
limit of 1Z < 25.5, we obtained a probability of <~ 0.95 
that the equivalent widths of the LBG- and narrowband- 
selected LAE samples could be described as originating 
from the same population. This result, consistent with 
Vcrhammc et al. (2008), is indicative of the similarity 
in equivalent width of objects spanning comparable UV 
luminosities. Conversely, the probability that the equiv- 
alent widths of our sample of LBG-selected LAEs and 
the faint (R > 25.5) Gronwall et al. (2007) LAEs derive 
from the same parent population is 8 x 10~ 8 . Similarly, 
when the Gronwall et al. (2007) data are divided using 
the boundary R = 25.5, a K-S test yields a probability 
of 7.9 x 10~ 6 that the two equivalent width populations 
derive from the same parent distribution. 

In other words, while the equivalent width distribu- 
tions of LBG- and narrowband-selected LAEs are quite 
disparate when considered at face value, it appears that 
the bulk of this disparity can be attributed to the differ- 
ent luminosities probed by the respective samples. When 
narrowband-selected LAEs are restricted to objects more 
luminous than R — 25.5, a K-S test reveals a high prob- 
ability that the equivalent width distributions of strong 
Lya-emitting LBGs and bright LAEs derive from the 
same parent population. This similarity in equivalent 
width, when a standard continuum magnitude range is 
adopted, is consistent with the picture that bright LAEs 
represent the same population as Lya-emitting LBGs 
(e.g., Verhamme et al. 2008). 

5.3. Escape of Lya Photons 

As Lya is susceptible to a variety of complex radia- 
tion transport effects including resonant scattering (e.g., 
Neufeld 1990) and the velocity field of the ISM (e.g., Ver- 
hamme et al. 2008), the escape of Lya emission is a probe 
of the physical conditions of galaxies' interstellar media. 
In the simplest case of a single-phase ISM in which gas 
and dust are well mixed, Lya photons experience more 
attenuation than continuum photons due to the longer 
dust absorption path lengths expected for resonantly 
scattered radiation. To explain the large Lya equiva- 
lent widths observed in some apparently dusty LAEs, 
the presence of a multi-phase ISM - in which dust is seg- 
regated in neutral clouds - has been proposed (Neufeld 
1991; Hansen & Oh 2006; Finkelstein et al. 2009). In 
such a geometry, Lya photons are scattered at the sur- 
face of dusty clumps while continuum photons penetrate 
the clumps and are preferentially scattered or absorbed. 
Galaxies in which the dust is segregated in neutral clouds 
are consequently predicted to exhibit larger Lya equiv- 
alent widths than would be expected given their stellar 
populations. 

We considered the likelihood of the above dust geom- 
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Fig. 12. — Lya escape fraction as a function of dust extinction, 
for objects with best-fit stellar ages t+ > 100 Myr. A significant 
correlation is observed, in the sense that objects with lower dust 
attenuations have larger escape fractions. The best-fit line from 
Verhamme et al. (2008) is overplottcd as a dashed line. 

etry in our sample of LBGs. Our methodology relied 
on comparing observed and predicted Lya luminosities, 
where the former, £ Lya (fl ux ), was calculated simply from 
the Lya line flux and the latter, iLya(SFR)) was derived 
from the conversion between SFR (from population syn- 
thesis modeling) and Lya luminosity, assuming the con- 
ventions (Case B recombination; T c = 10,000 K) of Kcn- 
nicutt (1998) and Brocklehurst (1971): 



^Lya(SFR) 



SFR [M Q yr~ 
9.1 x 10- 43 



(2) 



where SFR was calculated assuming a Salpeter (1955) 
IMF. Two similar parameters incorporating observed and 
predicted Lya luminosities were of interest: the escape 
fraction, f esc , and the relative escape fraction, f eS c,rel- 
The escape fraction is the ratio of observed to predicted 
Lya luminosities, whereas the relative escape fraction is 
this same ratio with an extra term corresponding to a 
dust correction at the rest wavelength of Lya. f e sc,rel 
is a probe of the degree to which Lya and continuum 
photons experience the same level of dust attenuation. 
In the case that f eS c,rei — 1, Lya and continuum photons 
are attenuated by the same E(B-V). On the other hand, 
a relative escape fraction less than unity is indicative of a 
dust geometry in which Lya photons are attenuated more 
than continuum photons. Conversely, a galaxy in which 
dust is segregated in neutral clouds and where the Lya 
flux is consequently enhanced relative to the continuum 
flux would have f esc ,rel > 1- We present the equations 
for f esc and f esc ,rel below, where fc'(1216) parameterizes 
the Calzetti et al. (2000) starburst attenuation law at 
1216 A: 



fesc = : LyatttUXJ , (3) 

^Lya(SFR) 

/ eS c, re ^/e SC XlO°- 4E ( B - V ) fc '( 1216 ). (4) 

When calculating f esc and f eS c,rei, we limited our anal- 
ysis to only those objects with W^ya > A due to the 
unphysical nature of inferring a Lya luminosity from a 



negative Lya flux. We also required objects to have best- 
fit stellar ages older than 100 Myr in order to ensure the 
validity of Equation 2, which assumes continuous star 
formation for at least that duration (Kcnnicutt 1998). 
These criteria isolated a sample of 105 objects. We note 
that the questions raised about the nature of the dust 
extinction law for young objects (§4.2) are not of con- 
cern given this older sample. To correct for slit losses, 
we normalized the G magnitudes determined from pho- 
tometry and spectroscopy; the latter were calculated by 
multiplying the G filter transmission curve (4780/1100 
A) over the optical spectra. The Lya line fluxes were 
adjusted accordingly, where the median flux correction 
was an increase by a factor of <~ 1.7 (0.6 magnitudes). 
This procedure relies on the assumption that Lya emis- 
sion has the same spatial extent as broadband 940-1470 
A emission. While some recent work (e.g., Steidel et al. 
2000; Matsuda et al. 2004; Hayashino et al. 2004) has 
suggested that Lya and UV continuum emission may be 
spatially decoupled with Lya more extended than con- 
tinuum emission, in the absence of simultaneous high- 
resolution Lya and UV continuum imaging for our en- 
tire sample we assume for simplicity that the line and 
continuum emission are coincident. 

The calculated values of f esc vary from 0.00-0.85. In 
Figure 12, we plot E(B-V) versus f esc ; we find a strong 
correlation in the sense that objects with large dust at- 
tenuations tend to have small Lya escape fractions (rsR = 
-0.712, P S R = 0.0000). We note that this trend is con- 
sistent with the existence of LBGs in our sample exhibit- 
ing both Lya in absorption and red UV continua, where 
these particular objects were not included in this anal- 
ysis due to their negative Lya fluxes. Other authors 
have derived similar results, although with significantly 
smaller sample sizes (e.g., Verhamme et al. 2008; Atek 
et al. 2009). The scatter present in our larger data set - 
some of which may be due to systematics in estimating 
f esc - is a useful probe of the nature of Lya radiative 
transfer in LBGs. We further discuss a physical picture 
of LBGs consistent with these f esc results in §6. 

In order to investigate the nature of the gas and dust 
distributions for the objects in our sample, we derived 
the relative escape fraction of the 105 LBGs with W^ya 
> A and i* > 100 Myr. We found that 103 objects had 
fesc,rel < 1 and the sample as a whole was character- 
ized by (/ esc (rel)) = 0.27 (Figure 13). For comparison, 
Gronwall et al. (2007) and Nilsson et al. (2009b) report 
typical fesc,rei values of ~ 0.30 and ~ 0.60 in LAE sam- 
ples at z ~ 3.1 and z ~ 2.25, respectively. Even if we 
conservatively assume Lya slit losses for our sample ap- 
proaching a factor of two (Hayashino et al. 2004), these 
low relative escape fractions, well below unity, do not 
support the hypothesis in which dust is primarily seg- 
regated to neutral clouds in Lya-emitting LBGs. We 
note that the absorption of Lya photons by the IGM 
is not a plausible explanation for these low relative es- 
cape fractions; at z ~ 3, Lya line fluxes are predicted 
to experience an IGM-induced decrement of only ~ 20% 
(Madau 1995; Shapley et al. 2006). Furthermore, the ob- 
served symmetry of the Lya lines in our sample (Figures 
1 and 7) are indicative of minimal IGM absorption. Ad- 
ditionally, if one adopts a sub-solar metallicity (in con- 
trast to Zq used in this paper), the calculated values of 
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Fig. 13. — Histogram of f esc ,rel values, for objects with > 
100 Myr. Note that the majority of objects (103/105) have / esc re j 
< 1, where the sample as a whole is characterized by (f esc . re i) = 
0.27. These low relative escape fractions are indicative of Lya pho- 
tons experiencing more attenuation than non-resonance continuum 
photons. These results are not consistent with the hypothesis of a 
clumpy interstellar medium, as has been suggested to exist in some 
LAEs (e.g., Finkelstein et al. 2009). 

fesc,rei will be lower than those presented here due to 
the fact that the ratio of the ionizing photon luminosity 
to the observed UV luminosity increases with decreasing 
metallicity (e.g., Leitherer et al. 1999). This result of 
relative escape fractions below unity indicates that Lya 
appears to experience more attenuation than continuum 
photons, consistent with the physical picture of a homo- 
geneous distribution of gas and dust in which resonantly 
scattered Lya photons have longer dust absorption path 
lengths than continuum photons. Such a simple picture 
is likely not an accurate representation of the ISM of z 
~ 3 LBGs; indeed, the ISM of our own Milky Way is 
known to be highly inhomogencous (Dickey & Garwood 
1989). However, our results suggest a scenario in which 
dust and gas are well mixed among the different phases 
of an inhomogeneous ISM. 

6. SUMMARY AND CONCLUSIONS 

We have analyzed a sample of 321 optically-selected z 
~ 3 Lyman break galaxies in order to investigate the rela- 
tionship between stellar populations and Lya emission. 
The equivalent width of the Lya feature was robustly 
estimated from rest-frame UV spectroscopy and broad- 
band G1ZJK S + Spitzer IRAC photometry was used to 
conduct stellar population modeling to derive the key 
properties of age, extinction, star formation rate, and 
stellar mass for all objects with photometric coverage in 
at least one near- or mid-infrared passband (<~ 80% of the 
sample) . The limited luminosity range of LBGs enabled 
a controlled investigation of the nature of Lya emission 
within a sample of objects spanning similar luminosities. 
Given the luminosity-dependent trends in Lya equivalent 
width, stellar populations, and clustering (e.g., Gronwall 
et al. 2007; Lai et al. 2008; Adelberger et al. 2005b), this 
controlled study represents an improvement over previ- 
ous work. We used a variety of statistical tests to an- 



alyze the correlation between Lya emission and stellar 
populations in LBGs from the standpoint of comparing 
spectroscopic observations, photometric data, and best- 
fit stellar population parameters. The relative Lya es- 
cape fraction was used to probe the relative distributions 
of gas and dust in the objects' interstellar media. Below, 
we summarize our results. 

• Lya equivalent width is correlated with age, SFR, 
and E(B-V), respectively, in the sense that galaxies 
with strong Lya emission are older, more quiescent, 
and less dusty than their counterparts with weak or 
absent Lya emission. Taking into account the un- 
certainties on the best-fit population parameters, 
the probability of a null hypothesis (i.e., no cor- 
relation) was excluded at the 3-5ct level for these 
relationships. We found that stellar mass was not 
significantly correlated with W^ya ~ the null hy- 
pothesis could be excluded at only the <~ la level. 
These results were consistently derived from several 
different analyses of the data, including correlation 
tests, composite spectra, and the binary division 
of the sample into LAEs (WL ya > 20 A) and non- 
LAEs (W hya < 20 A). 

• Analysis of the relative escape fraction of Lya is 
consistent with Lya photons experiencing more 
attenuation than non-resonance continuum pho- 
tons. We also find that the Lya escape fraction 
is strongly correlated with E(B-V), where galaxies 
with more dust attenuation also have lower escape 
fractions. 

The observed correlations between Lya emission and 
stellar populations are consistent with the physical pic- 
ture proposed by Shapley et al. (2001), in which young, 
dusty LBGs experience vigorous outflows from super- 
novae and massive star winds. Shapley et al. (2003) re- 
ported evidence that both gas and dust are entrained in 
the outflows, evidence that these "superwinds" could ten- 
ably decrease a galaxy's dust and gas covering fraction 
over several tens of Myr. While the more mature LBGs 
may have as much or more overall dust content than 
younger galaxies, the fact that these older objects (i* > 
100 Myr) are also more likely to exhibit Lya emission is 
an important point: this trend demonstrates that it is 
the dust covering fraction - not the total amount of dust 
- that modulates Lya emission. Independent evidence 
for a lower covering fraction of dusty gas in older objects 
is provided by a trend toward weaker low-ionization in- 
terstellar absorption lines with increasing age in the com- 
posite spectra presented in §4.2. 

While this physical picture has been put forth previ- 
ously, our current results are more strongly supported 
due to the wide range of statistical tests and analysis 
methods that we employed. We compared our larger 
sample of data in different forms - spectroscopic, pho- 
tometric, and best-fit parameters - using a battery of 
statistical methods including correlation and K-S tests. 
Lya emission strength was furthermore treated as both 
an independent and dependent variable. These analyses 
represent a systematic, powerful approach to elucidating 
trends. 

Results of both the escape fraction and the relative 
escape fraction can be used to further constrain the 
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physical picture of LBGs and bright narrowband-selected 
LAEs. While the observed inverse relationship between 
E(B-V) and f esc is consistent with earlier results (e.g., 
Verhamme et al. 2008; Atck ct al. 2009), the scatter 
present in our larger sample highlights the myriad factors 
modulating the Lya escape fraction (e.g., galaxy kine- 
matics, dust, and outflow geometry). Further analysis of 
this scatter will be instrumental in constraining how the 
Lya escape fraction varies. The observed relative escape 
fraction of Lya emission - where the majority of LBGs 
have fesc,rei below unity - is indicative of a physical pic- 
ture in which LBGs contain gas and dust that is well 
mixed. In this geometry, Lya photons experience more 
attenuation than continuum photons do because of the 
former's increased path lengths from resonant scattering. 
This is in contrast to the scenario proposed by Finkel- 
stein et al. (2009), in which dust is segregated in neutral 
gas clumps surrounded by an ionized, dust-free medium. 
We caution, however, that the relative escape fraction of 
Lya is likely modified by a variety of factors including 
galaxy kinematics, dust, and outflow geometry. As such, 
additional observations are necessary in order to more 
fully understand the distribution (and evolution) of gas 
and dust in LBGs. 

While we hope that this work has illuminated the rela- 
tionship between Lya emission and stellar populations in 
2^3 LBGs, future systematic studies at other redshifts 
- ideally over larger dynamic ranges in UV luminosity - 
are necessary in order to address the discrepancies be- 
tween trends reported here and those at z <~ 2 or z <~ 
4 (e.g., Erb et al. 2006a; Pentericci et al. 2007). The 
nature of complex Lya emission morphologies - such as 
the double-peaked objects that comprise ~ 4% of the 
sample discussed here - has also yet to be explored in 
light of understanding if these objects' stellar popula- 
tions and gas distributions differ from those of more typ- 
ical Lya-emitters or absorbers. Similarly, an analysis of 
the geometry of dust and gas in the interstellar media 
of z <~ 3 LBGs is vital to understanding what kind of 
dust attenuation law - the Calzetti et al. (2000) relation, 



a SMC-like law, or another relation completely - most 
accurately describes these objects. Furthermore, a basic 
property critical to placing galaxies in an evolutionary 
context - metallicity - remains poorly constrained for 
the z <~ 3 LBG population (but see Maiolino et al. 2008). 
While metallicity has been investigated in the rare cases 
of gravitationally lensed objects (e.g., Teplitz et al. 2000; 
Pettini et al. 2002; Hainline et al. 2009; Quider et al. 
2009), a systematic study of elemental abundances in a 
larger sample of z ~ 3 LBGs would be illuminating (akin 
to the z ~ 2 sample discussed in Erb et al. 2006a). Given 
anticipated advances in both telescopes and instrumen- 
tation, we hope that these questions are illustrative of 
many others that will help to guide future LBG research. 
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